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I.  INTRODUCTION 


Tany  scientific  and  technical  endeavors  require  the  reconstruct  ion 
of  a  three-dimensional  solid  from  a  collection  of  two-dimensional 
planar  contours.  These  contours  are  obtained  by  some  sensor  method 
that  samples  the  original  three-dimensional  solid  along  a  finite  number 
of  parallel  planes.  The  data  extracted  from  that  set  of  parallel 
planes  are  contours  that  lie  along  the  solid’s  exterior  and  interior 
surfaces.  ^he  contours  on  the  parallel  planes  appear  as  line  segments, 
■^he  line  segments  are  either  closed  loops,  open  segments,  or  single 
points.  The  goal  of  surface  construction  is  the  formation  of  surface 
patches  between  contours  on  adjacent  planes  such  that  an  approximation 
of  the  original  three-dimensional  solid  is  formed. 

Surface  construction  fay  the  triangulat ion  of  two-dimensional 
contours  is  the  procedure  by  which  a  pair  of  parallel,  planar  contours 
are  "mapped  together"  and  then  "triangulated"  into  surface  patches  that 
form  a  surface  display.  The  mappirig  operation  of  the  surface 
construction  algorithm  identifies  which  contours  on  consecutive, 
parallel  planes  should  be  mapped  together,  and  exactly  which  portions 
■;f  those  contours  should  be  connected.  The  triangulat ion  operation 
“^orms  the  connections  between  contours  'On  adjacent  planes  by  building 
triangular  tiles  between  those  mapped  contours.  Each  triangular  tile 
IS  b  r.lt  from  an  .ndividual  line  segment  from  one  contour  and  a  single 
point  '^ri'M  the  end  of  a  line  segment  on  the  other  mapped  contour.  This 
t.ling  -operation  is  performed  for  all  line  segments  in  the  connect 


region  of  each  mapped  contour.  The  connect  region  is  that  section 
of  coordinates  designated  as  mappable  for  a  pair  of  contours  on 
consecutive  planes. 

N'otat  ional  ly,  this  problem  has  been  specified  as  follows: 

"fin  unknown  three  dimensional  solid  is  intersected  by  a  finite 
number  of  specified  parallel  planes. 


only  information  about  the  solid  consists  of  the 
; ntersect ions  of  its  surface  with  the  planes.  Each  of  these 
•. ntersect  ions  is  assumed  to  be  a  simple  closed  curve.  These  curves 
are  not  completely  specified;  instead,  a  finite  sequence  of  points 
encountered  during  a  positive  (counterclockwise)  traversal  of  each  of 
the  original  curves  is  given.  The  curve  segment  between  two 
consecutive  points  is  approximated  by  a  linear  segment,  called  a 
contour  segment. 


VJe  reduce  the  problem  of  constructing  such  an  approximating 
surface  to  one  of  constructing  a  sequence  of  partial  approximations, 
each  of  them  connecting  two  contours  lying  on  consecutive  planes. 
[Figure  l.i: 

_et  one  contour  be  defined  by  the  sequence  of  rn  distinct  contour 
points  P2i,  PI,  ...,  P(m-l),  and  let  the  other  contour  be  defined  ay 
the  sequence  of  n  distinct  contour  points  D0,  Ql,  ...,  Q(n-l).  We 
note  that  P0  follows  P(m-l)  and  that  Q0  follows  Q(n-l),  and  so 
indicies  of  P  are  modulo  m  and  indicies  of  Q  are  modulo  n.  We  wish 
to  create  a  surface  between  the  contours  P  and  Q.  The  surface  is 
constructed  of  triangular  tiles  between  these  two  contours.  The 
verticies  of  these  tiles  are  contour  points,  with  the  verticies  of 
each  tile  tal/en  two  from  one  sequence  and  one  from  the  other.  Thus, 
each  tile  is  de'^ined  by  a  set  of  three  distinct  elements  either  of 
the  fo^'m  •CPi,Pk,Cj/  or  ■Cai,Qk,Pj>,  [Figure  1.23 

Each  tile’s  boundary  will  consist  of  a  single  contour  segment  and 
two  scans,  each  connecting  an  end  of  the  contour  segment  with  a 
i:vmon  point  on  the  other  contour. "  [Ref,  13 

“"his  nctational  specif icat ion  of  the  problem  is  consistent  in  all 
pacers  accessible  in  the  literature  on  surface  construction  CRef.  13 


["ef. 


23  CRef.  33  CRef.  43. 


The  initial  emphasis  of  this  paper  is  a  review  of  the  previous 
algorithms  for  surface  construction.  Included  in  this  review  is  a 
discussion  of  each  algorithm’s  capabilities  and  limitations.  After 
this  review,  we  present  a  new  algorithm  for  surface  construction  that 
is  more  comprehensive  than  any  that  has  previously  appeared  in  the 
literature.  -ollowing  that  discussion,  we  examine  the  limitations  of 


our  new  algorithm. 


II.  LITERATURE  REVIEU 


Ir>  order  to  understand  the  nature  of  our  new  algorithm  for  surface 
construction,  it  is  irnoortant  that  we  understand  the  previous 
algorithms  for  surface  construction.  Four  such  algorithms  have 
provided  the  background  necessary  for  the  development  of  our  algorithm. 
These  are  presented  in  chronological  order. 

A.  FUCHS  ALGORITHM 

'I'he  first  algorithm  we  examine  for  the  reconstruct  ion  of  a 
three-dimensional  object  from  its  planar  contours  was  presented  by 
Fuchs  in  1977  [Ref.  11.  The  problem  statement  from  that  article 
(reproduced  in  our  introduction)  has  been  used  in  all  subsequent  papers 
which  build  upon  the  Fuchs  algorithm.  The  major  contribution  of  that 
article,  in  addition  to  the  concise  statement  of  the  problem,  is  the 
presentation  of  an  algorithm  capable  of  connecting  simple,  closed 
contours  (Figure  2.1). 

The  problem  with  the  Fuchs  algorithm  stems  from  its  inability  to 


hand  le 

multiple 

C'ontours  on 

adjacent  planes 

(Figure 

2.2). 

i  1 1  onal  1  y ,  rio 

mechanism  is 

provided  to  handle 

partial 

contour 

tiiap'iiings 

■:*r  open 

(non-closed ) 

contours.  With  respect  to  the 

case  of 

■  contours  on  adjacent  planes,  no  mechanism  is  provided  to 
..entify  which  of  the  contours  should  be  mapped  together.  The  general 
casa  for  surface  construction  is  to  have  multiple  contours  on  each 
plane.  The  problem  with  partial  contour  mappings  is  that  the  Fuchs 
alj'-'  ithn  can  .only  construct  a  complete  triangulation  between  adjacent 


The  Christiansen  algorithm  is  not  capable  of  handling  open  contours, 
nor  is  it  capable  of  handling  complex  cases  of  multiple  contours  on 
adjacent  planes,  except  by  way  of  expensive  human  interaction.  ft  final 
note  of  interest  with  respect  to  this  algorithm,  is  the  use  of  a 
heuristic  for  selection  of  the  nodal  connections.  In  cases  where 
contours  on  adjacent  planes  are  mutually  centered  and  are  reasonably 
similar  in  size  and  shape,  selection  for  nodal  connection  is  based  on 
"shortest  diagonal"  rather  than  minimum  triangular  area  CRef.  £:  p. 
1381.  During  this  operation,  one  of  two  nodes  is  selected  to  create 
the  next  triangular  surface  patch.  The  nodes  under  considerat ion  are 
the  two  "next"  nodes  of  each  contour.  By  determining  the  length  of 
each  of  the  possible  diagonals  for  the  surface  patch,  the  connection 
'■ode  is  selected  based  on  minimum  length. 

C.  SnftNTZ  ftLSQRlTHV, 

"he  algorithm  presented  in  the  Shantr  article  [Ref.  31  extends  the 
algorithms  of  Fuchs  and  Christiansen  to  handle  contour  defined  objects 
vh.ch  are  highly  branched  and  have  holes.  Multiple  contours  on 
adjacent  planes  are  handled  by 

.  .  first  concatenating  the  contours  on  each  plane  into  a 

single  large  contour  using  minimum  distance  links,  then  performing 
the  mapping  between  the  resulting  composite  contours."  CRef.  3:  d. 
2421 

Fhantz  uses  the  simple,  closed  contour  mechanism  of  Fucns  to  form  tne 
connections  between  the  composite  contours.  Once  the  connections  have 
been  formed,  the  extraneous  ones  (due  to  concatenat ion)  are  removed. 
Some  difficult  multiple  contour  cases  for  this  algorithm  require  human 


inter'action  to  solve  ambiguities.  Similar  to  the  Christiansen 
algorithm,  Shantz  states  that  this  should  be  avoided  since  human 
interaction  is  "extremely  labor  intensive."  He  cites  a  case  which 
required  50  to  30  hours  of  contour  splitting,  using  an  interactive 
cursor,  to  produce  a  surface  display  for  the  highly  convoluted  cortex 
and  basal  ganglia  contours  (extracted  from  the  Livingston  brain 
catabase) • 

This  algorithm,  similar  to  the  Christiansen  algorithm,  is  limited 
in  its  ability  to  handle  cases  of  open  contours  and  partial  contour 
mappings.  Also,  cases  of  multiple  contours  on  adjacent  planes  can  be 
handled  only  when  a  composite  contour  can  be  formed,  or  when 
ambiguities  are  resolved  via  human  interaction. 

GPMfiPOTHY  ALECRITHM 

The  most  recent  algorithm  for  surface  construction  from  planar 
contours  was  presented  in  a  paper  by  Ganapathy  CRef,  4].  That 
algorithm  is  essentially  an  improvement  on  the  Fuchs  and  Christiansen 
algorithms  for  simple,  closed  contours,  without  the  capabilities 
tesc-'ioed  by  Shantz.  Like  Fuchs,  Ganapathy  assumes  a  complete  mapping 
of  contours,  which  is  not  always  possible.  The  improvement  over  the 
-uchs  and  Ch/r i st i ansen  algorithms  is  attributed  to  the  use  of  a  more 
computat ional ly  expedient  heuristic  for  triangulations. 

The  problem  with  the  Ganapathy  algorithm  is  that  it  presents  a 
general  solution  for  handling  only  the  simple  case  of  mapping  single, 
closed  contours  on  adjacent  planes.  The  issues  of  multiple  contour 
mappings  and  partial  contour  mappings  are  ignored.  Additionally,  no 


mechanism  for  user  interaction  is  provided  for  resolving  mapping 
ambiguities,  further  limiting  the  algorithm  to  simple  cases. 

E.  SuMMRRY 

!Mone  of  the  above  papers  provides  a  complete  solution  to  the 
problem  of  surface  construction  via  the  triangulation  of  contours, 
'.x'hat  is  required  is  an  algorithm  with  capabilities  for  multiole 
contours  per  plane  and  partial  contour  mappings.  Additionally,  the 
algorithm  should  suoport  simple  cases  of  branching  and  provide  a 
mechanism  for  human  interaction  for  the  resolution  of  highly  ambiguous 
mappings. 

The  surface  construction  algorithm  we  present  handles  not  only  the 
simole  contour  mapping  problem,  but  also  provides  a  more  comprehensive 
procedure  for  solving  the  multiple  contours  per  plane  and  partial 
mapping  problems.  The  only  capability  lacking  from  our  algorithm  is 
that  for  handling  branching  as  per  the  Christiansen  paper.  A  detailed 
discussion  of  our  algorithm  follows,  with  a  presentation  of  algorithm 


■euristics  and  limitations. 


T  ▼  T 
i  i  A  • 


THE  ALGORITHM 


Iio  the  last  section,  we  presented  a  discussion  of  previous 
algorithms  for  surface  construction  via  triangulat ion.  Here,  we 
present  an  in-depth  discussion  of  our  algorithm  by  first  discussing 
'■'nown  inout/output  data  structures.  Following  this  presentation,  an 
overview  of  the  major  parts  of  the  algorithm  precedes  a  detailed 
discussion  of  the  parts. 

a.  INPUT/QUTPUT  SPECIFICfiTIONS 

"he  problem  of  surface  construction  of  an  object  from  a  set  of 
planar  contours,  as  seen  in  Figure  3.1,  can  be  reduced  to  one  of 
constructing  the  surface  tr langulat ions  between  two  adjacent  planes, 
"he  speci f icat ion  of  the  problem  can  be  best  viewed  by  detailing  the 
-(.'cwn  input  data  structures: 

*  total (i)  :  number  of  contours  on  plane  i. 

^  start(j,i)  :  start  of  contour  j  on  plane  i. 

*  length(j,i)  ;  number  of  coordinates  in  contour  j  on  plane  i. 

*  type(j,i)  :  tyoe  of  contour  j  on  plane  i. 

(CLaSED_LOaP,  □PEN_SEG!»‘ENT,  or  SINGLE_PDINT) 

*  interior (j, i )  :  value  of  contour  j’ s  interior  with  respect  to 

the  contour  line. 

(HIGH,  LOW,  or  INDETERM  I NOTE) 

*  coords ( XYZ, pointer, 1 ) :  input  coordinates  for  all  contours  on  plane 

1.  To  isolate  contour  j  on  plane  i: 

for  (pointer  =  start (j,i)  +  k  -  1), 

where  k  =  1,  length(j,i). 

-rom  the  above  data,  we  desire  to  produce  the  following  output  data 
struct  ures : 


*  nuni_coords  ;  number  of  coordinates  generated  for  the  two 

input  planes. 

*  new_coords (XYZ, num_coords) :  coordinates  generated  by  the  surface 

construction  process  for  the  two'  planes. 

*  new_conns (nufn_coord5) :  drawing  instructions  for  each  coordinate 

generated  (SETPOINT,  DRfiWTO,  DRfiWPOINT). 

If  the  output  data  is  in  the  form  of  triangular  surface  patches,  an 
alternative  data  structure  is  required; 

*  nurn_patche5  :  number  of  surface  patches  generated  for  the 

input  two  planes. 

*  new_coords ( XYZ )  :  new  coordinates  generated  by  the  connection 

process. 

*  patches (3, num_patches) :  a  3  by  num_patches  array  of  triangles. 


B.  THE  fiLGCRITHM 

Our  surface  construction  algorithm  is  composed  of  the  following 
outlined  steps: 

; 1 )  Input  and  Inventory  Compilation; 

The  data  structures  defining  the  contours  are  processed  to 
extract  the  pertinent  data.  This  data  includes  the  number  of 
contours  per  plane,  the  coordinates  defining  those  contours  and 
the  types  of  the  contours.  Pddit ional ly,  two-dimensional 
bounding  boxes  are  described  about  each  contour  for  processing 
considerat ion  in  step  2.  This  compilation  of  data  creates  the 
data  structures  required  for  surface  construction. 

(3)  Overlap  Determination  and  Contour  Item  tapping; 

In  this  step  of  the  alg-orithm,  we  determine  which  contours  on 
adjacent  planes  have  significant  overlap,  and  which  contoui's’ 
exteriors  are  near.  This  information  is  used  to  designate  which 


contours  should  be  connected  via  tniangulat ions.  The  assignment 
of  overlap  is  accomplished  through  the  use  of  a  value  for  the 
overlap  percentage.  This  value  is  computed  from  the  areas  of 
the  two-dimensional  bounding  boxes,  as  seen  in  Figure  3.2,  of 
each  contour.  The  overlap  percentage  is  used  to  give  priority 
to  contour  mappings  that  have  the  highest  percentage  of  total 
overlap  area. 

In  this  step  of  the  algorithm,  we  also  perform  consistency 
checks  for  each  contour  pair.  One  such  consistency  check  is 
executed  using  the  contour  interior  specification  and  the 
overlap  percentage  value.  Contour  interior  specifications  are 
assigned  as  the  value  of  a  contour  with  respect  to  its  immediate 
interior.  fls  such,  a  contour  is  LOW  valued  if  it  is  taken  from 
the  exterior  of  a  solid  object,  such  as  the  skin  of  an  apple. 
Conversely,  a  contour  is  HIGH  valued  if  its  immediate  interior 
is  non-solid.  Using  these  pieces  of  information,  we  are  able  to 
eliminate  contour  mappings  of  high  overlap  percentage  which 
'esult  in  erroneous  approximations  of  the  original 
three-dimensional  solid. 

illustrate  the  application  of  this  consistency  check,  let  us 
consider  the  mapping  example  for  Figure  3.3.  Here  we  are 
presented  with  a  set  of  contours  taken  fi'om  a  solid  cone 
standing  within  a  hollow  cone.  In  this  case,  contour  1  on  plane 
1  has  a  high  overlap  percentage  with  contour  2  on  plane  2. 
However,  since  contour  2  on  plane  2  is  low  valued  with  respect 


contour  1 ,  plane  2 


to  its  solid  iriterior  and  contour  1  on  plane  1  is  high  valued, 


this  mapping  can  be  eliminated. 

The  interior  specifications  are  also  used  to  determine  whether 
the  mapping  is  interior  to  interior  or  exterior  to  exterior.  fln 
interior  to  interior  mapping  is  one  which  maps  the 
interior  of  one  contour  to  the  interior  of  another  contour, 
■'his  form  of  mapping  is  indicative  of  contours  taken  from  a 
surface  with  a  shallow  gradient,  i.  e. ,  a  surface  where  the 
mapped  contours  are  of  similar  size  arid  shape,  and  where  the 
contours  have  significant  overlap.  An  exterior  to 
exter lor  mapoing  is  one  which  maps  the  exterior  of  one 
contour  to  the  exterior  of  another  contour.  This  form  of 
maoping  is  indicative  of  contours  taken  from  a  surface  with  a 
steep  gradient,  i.  e. ,  a  surface  where  mapped  contours  are  of 


dissirni  lar 

size  and 

shape,  and 

where 

the  contours 

overlap 

percentage 

IS  slight. 

Interior 

to  interior  mappings 

are  more 

oO'fiimoin.  T 

he  exterior 

to  exterior 

mapping 

is  indicated 

for  cases 

of  two  contours  with  a  low  percentage  of  overlap  and  differing 
Interior  specifications  (HIGH; LOW,  or  vice  versa). 

-’irv  the  Coordinate  irapDino  for  each  flapped  Contour  Pair: 

r  each  coordinate  pair  fro'm  step  two,  we  form  a  complete 
ordinate  to  coordinate  mapping.  A  coordinate  mappino 

a  tentative  set  of  triangulat ion  connections  oetween  the 
rtour  pairs.  There  are  two  procedures  for  determining  this 
it.al  coordinate  mapoing.  The  procedure  used  is  dependent  on 
a  type  ijf  mapoing  found  for  the  paired  contours  in  the 


previous  step  (interior  to  interior,  or  exterior  to  exterior). 


Oddit ional ly,  both  procedures  try  to  form  triangulation  segments 


of  shortest  length,  similar  to  the  Christiansen  algorithm.  fi 


general  statement  of  this  selection  process  is  that  we  are 


trying  to  map  coordinate  i  of  contour  n,  plane  1  to  coordinate  j 


of  contour  rn,  plane  £  such  that  the  distance  between  the  two 


coordinates  is  minimised.  On  additional  qualification  toi  this 


distance  minimising  criterion  is  that  coordinate  connections  go 


not  cross,  i.  e. ,  coordinates  3  and  4  of  plane  1  are  not  maoped 


to  coordinates  6  and  5  of  plane  £  respectively. 


(4)  Cont inuit V  Recocnit ion; 


The  coordinate  to  coordinate  mapping  formed  in  step  three  is 


examined  for  continuity.  Continuity,  in  this  case,  is  defined  as 


I'llows.  -irst,  we  form  sets  of  coordinates  from  the  coordinate 


mapping  such  that  each  coordinate  of  each  set  is  constrained 


within  a  coordinate  tolerance  and  within  a  distance  range.  The 


coordinate  tolerance  factor  is  a  ratio  of  the  number 


5f  coordinates  in  the  larger  contour  divided  by  the  number  of 


coordinates  in  the  smaller  contour  times  a  window  value.  The 


tolerance  factor  is  used  to  group  coordinates  into  a  single  set 


based  upo'n  their  mapped  coordinate  number  being  within  plus  or 


minus  tolerance  of  the  last  mapped  coordinate  added  to  the  set. 


'he  tolerance  sets  formed  are  then  compared  for  overlapping 


distance  ranges.  Any  sets  that  have  overlapping  distance  ranges 


are  then  merged.  The  merged  set  with  the  smallest  distance  in 


4 


I  I 
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it  is  the  set  of  coordinates  for  which  connections  should  be 


generated,  fill  other  coordinates  are  left  unconnected. 

(5)  i^apaino  Cancellation; 

Cnee  we  have  decided  to  generate  the  connections  for  a  part  of  a 
contour,  we  cancel  any  further  mappings  to  that  piece  of  the 
contour.  This  operation  is  required  for  partial  mappings  in 
which  two  or  more  contours  on  one  plane  are  to  be  mapped  to  a 
single  contour  on  another  plane.  This  cancellation  precludes 
connecting  contour  points  which  have  already  been  selected  for 
connect  ion. 

(B)  Coinnection  Formation; 

Ue  generate  the  coordinates  for  the  triangulation  connections 
specified  in  step  four.  "In  between"  coordinates,  coordinates 
not  directly  mapped  but  within  the  tolerance  factor  for  the 
connection  mapping,  are  also  added  to  the  picture.  The  goal  of 
the  process  is  to  form  minimum  area  triangular  surface  patches 
for  each  segment  of  the  mapped  connection  region. 

1 .  Input  and  Inventory  Coimpi  lat  ion 

The  input  data  to  the  algorithm  consists  of  the  contour 
descriptions  '^or  two  adjacent  planes  of  a  three-dimensional  solid.  The 
purpose  of  this  step  of  the  algorithm  is  to  segment  this  data  into 
separate  contour  descriptions  and  to  determine  the  individual 
characteristics  of  each  contour.  Figure  3.4  consists  of  two  adjacent 
planes,  each  having  three  concentric  rings  of  similar  shape  and 
continuity.  Figure  3.6  consists  of  two  closed  loops  on  each  of  its 


Fig.  3.6  -  Example  of  a  set  of  contours  requiring  partial  mappings 
and  an  exterior  to  exterior  mapping;  (1,1)  and  (2,1)  to  (2,2). 


-  HIGH  interior  value 
LOW  interior  value 


Fig.  3.7  -  Connection  of  Figure  3.6,  with  contour  interior  values 
for  each  contour. 


planes. 


lane  1  has  two  small  interior  lobes,  while  plane  2  has  one 
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lar'ge  surrounding  contour  with  a  small  interior  contour.  The  contour 
descriptions  for  these  figures  are  composed  of: 
the  starting  coordinate  location, 
the  total  number  of  coordinates, 
the  contour  types, 
the  interior  values,  and 

the  contours’  two-dimensional  bounding  boxes. 

With  the  exception  of  the  interior  values,  all  of  these  characteristics 
are  easily  obtainable  from  the  input  data. 

The  procedure  necessary  to  obtain  the  contour  interior 
specifications  requires  an  evaluation  of  the  data  values  lying  along 
and  interior  to  the  contour  (see  Figure  3.3).  If  these  values  are  not 
contained  in  the  input  data,  a  mechanism  needs  to  be  provided  to  allow 
for  user  specification  of  contour  interior  values.  The  range  of 
interior  values  is  HIGH,  LOW  or  INDETERMINfiTE.  Without  this  value  the 
contour  pairing  operation  encountered  in  the  multiple  contours  per 
plane  situation  is  difficult.  In  that  case,  some  form  of  human 
interaction  is  necessary  to  designate  which  pairs  of  contours  should  be 
mapped  together.  If  an  interior  value  is  not  available,  and  the 
mapping  situation  is  not  complex,  it  can  be  set  to  INDETERMINRTE 
without  surface  construction  degradation. 


2.  Overlap  Determination  and  Contour  WaoDinn 

The  overlap  determination  and  contour  mapping  procedure  of  the 
surface  construction  algorithm  is  the  process  by  which  tentative 


contour  to  contour  mapping  assignments  are  made.  The  contour 
characteristics  which  are  necessary  for  this  procedure  are  the 
two-dimensional  bounding  boxes  and  the  contour  interior  specifications. 
This  mapping  process  is  the  key  component  in  the  disambiguation  of 
multiply  paired  contours. 

The  overlap  determinat ion  and  contour  mapping  procedure  is 
accomplished  in  the  following  manner.  First,  the  two-dimensional 
bounding  box  of  each  contour  on  plane  1  is  compared  for  overlap  with 
the  two-d imensional  bounding  box  of  each  contour  on  plane  2.  The 
coordinates  which  define  these  bounding  boxes  are  the  minimum  and 
maximum  X  and  Y  coordinates  from  each  of  the  contour  descriptions. 
(Additionally,  these  coordinates  are  adjusted  by  a  constant  value  to 
promote  overlap  for  exterior  to  exterior  mapping  situations.)  From 
this  operation,  a  table  called  the  overlap  table  is  produced.  It  is  a 
two-dimensional  table  that  contains  a  value  for  each  possible  pairing 
of  contours  between  the  two  planes.  The  value  recorded  in  each  table 
entry  indicates  the  extent  to  which  each  contour  overlaps.  If  there  is 
ro  bounding  box  overlap  for  a  pair  of  contours,  a  value  of  0.0  is 
recorded  in  the  table.  If  there  is  overlap,  the  value  recorded  in  the 
‘able  represents  the  percentage  of  overlap  with  the  larger  of  the  two 
contours.  "^his  value  is  computed  by  dividing  the  area  of  the  bounding 
box  overlap  by  the  area  of  the  bounding  box  of  the  larger  contour. 

After  the  overlap  percentage  has  been  computed  for  a  contour 
pairing,  it  is  used  in  conjunction  with  the  interior  specifications  to 
determine  the  mapping  type  for  the  contour  pair.  An  interior  to 
interior  mapping  is  indicated  when  a  high  percentage  of  overlap 


(greater  than  10%)  exists  for  a  pair  of  contours,  fl  consistency  check 
for  matching  interior  specifications  is  performed  for  every  pair  of 
contours  that  exhibits  this  high  an  overlap.  The  consistency  check 
requires  that  each  contour  pair  have  either  HIGHrHIGH,  L0W:L0U,  or 
INDETERMINPTE; anything  (HIGH  or  LOW)  interiors.  Contour  pairings  with 
high  overlap  but  inconsistent  interior  specifications  result  in  an 
adjustment  to  the  overlap  table  of  0.0  percentage  of  overlap.  fin 
exterior  to  exterior  mapping  is  indicated  when  the  overlap  percentage 
15  low  (less  than  10%)  and  item  interiors  are  non-matching.  Finally, 
all  contours  with  low  overlap  percentages  and  matching  interiors  are 
zeroed  in  the  overlap  table. 

Figures  3.3  and  3.9  graphically  represent  the  overlap 
determination  and  contour  mapping  for  Figures  3.4  and  3.6.  Included  in 
these  'igures  are  the  overlap  tables  produced  by  this  procedure.  The 
table  in  Figure  3.8  shows  three  valid  overlap  percentages  for  three 
different  contour  pairs!  (1,1)  -  (1,2),  (2,1)  -  (2,2),  and  (3,1)  - 

(3.2) .  '^our  of  the  entries  have  been  zeroed  by  the  consistency  check 
mechanism.  Without  this  capability,  high  valued  overlap  percentages 
wtould  appear  in  the  overlap  table  with  human  interaction  required  for 
their  disambiguation.  The  table  in  Figure  3.9  shows  two  high  overlap 
percentages  and  two  low  overlap  percentages.  This  data  indicates  that 
contours  (1,1)  and  (2,1)  both  map  interior  to  interior  with  contour 

(1.2) .  The  low  overlap  percentages  indicate  that  contours  (1,1)  and 
(2,1)  map  exterior  to  exterior  with  contour  (2,2). 
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The  coordinate  mapping  formation  procedure  for  each  coordinate 
pair  having  a  non-zero  overlap  (in  the  overlap  table)  begins  with  the 
pair  having  the  largest  overlap  percentage.  All  remaining  steps  in  the 
surface  construction  algorithm  are  carried  out  on  this  pair  before  the 
next  pair  of  contours  is  considered  for  mapping.  Mapping  paired 
contours  is  on  a  largest  to  smallest  overlap  percentage  ordering. 
Since  exterior  to  exterior  mapoings  are  indicated  only  in  situations 
where  the  overlap  percentage  is  low,  they  are  considered  for  mapping 
only  after  all  interior  to  interior  mappings  have  been  performed.  This 
study  follows  that  ordering  and  completes  the  description  of  the 
interior  to  interior  mapping  process  before  considering  the  separate 
process  necessary  for  exterior  to  exterior  mappings. 

“he  first  operation  performed  on  an  interior  to  interior 
overlap  pair  is  the  determination  of  which  contour  is  interior  to  the 
other.  This  assignment  is  accomplished  by  comparing  bounding  box  areas 
^or  the  contour  pair  and  designating  the  contour  with  the  smaller  area 
as  interior.  Once  the  interior  contour  assignment  has  been  made,  the 
center  coordinate  of  that  contour’s  bounding  box  is  computed. 

^'-e  knowledge  of  the  center  coordinate  of  the  interior  contour 
IS  used  in  the  following  manner.  For  each  coordinate  of  the  inner 
contour,  we  determine  which  coordinate  of  the  outer  contour  is  closest 
to  a  vector  drawn  from  the  center  coordinate  of  the  inner  cont>our 
through  the  coordinate  of  the  inner  contour  (see  Figure  3.10).  We  add 
the  qual i f icat ion  that  the  outer  coordinate  selected  by  this  procedure 
iiust  be  farther  ^rom  the  center  coordinate  than  the  inner  coordinate. 


Also,  the  outer  coordinate  must  be  on  the  same  side  of  the  vector  as 
the  inner  coordinate.  The  outer  coordinates  selected  by  this  mapping 

process  are  recorded  as  the  tentative  coordinate  mao 

coordinate  for  each  inner  coordinate.  We  also  record  the 
two-dimensional  distance  from  each  inner  coordinate  to  its  tentatively 
mapped  outer  coordinate.  The  resulting  data  structure  contains  the 
mapped  outer  coordinates  with  the  distance  to  the  inner  coordinate  to 
which  it  is  mapped. 

The  tentative  connection  map  for  Figure  3.4  is  very  good.  Due 
to  the  similarity  in  size  and  shape  of  the  mapped  contour  pairs,  there 
is  very  little  variation  in  the  mapped  distance  values  and  the 
coordinates  selected  for  mapping  appear  sequential.  On  the  other  hand, 
it  can  be  seen  in  Figure  3. 11,  that  large  variations  in  distance  values 
result  from  this  tentative  mapping  process,  and  mapped  outer 
coordinates  appear  with  large  gaps  in  the  sequencing.  This  is  due  to 
the  dissimilarity  of  the  contour  pair;  the  inner  contour  is  relatively 
simple  and  much  smaller  than  the  convoluted  outer  contour.  The 
procedure  used  to  delineate  a  correct  mapping  from  this  tentative 
mapping  is  described  below. 

a.  Continuity  Recognition 

the  continuity  recognition  procedure  uses  the  tentative 
connection  map  and  associated  distances  for  a  pair  of  contours  to 
ilsterrnine  the  set  of  coordinate  mappings  that  should  be  made  for  that 
pair.  In  the  previous  step  of  the  algorithm,  we  produced  the  tentative 
connection  map  for  all  of  the  coordinates  of  the  inner  contour.  This 
provides  a  rough  approximation  of  the  final  mapping,  but  it  must  be 


noted  that  not  all  of  the  inner  coordinates  need  be  involved  in  the 
final  mapping  for  that  pair.  The  continuity  recognition  procedure 
builds  sets  of  coordinate  mappings  that  are  both  continuous  and  of 
similar  mapped  distance  range.  These  continuity  sets  are  then  used  to 
determine  the  coordinate  sequences  that  should  comprise  the  final 
connection  mapping. 

The  first  step  in  this  procedure  is  to  assign  each 
coordinate  pairing  of  the  tentative  connection  map  to  an  initial 
continuity  set.  This  is  accomplished  by  stepping  through  the 
coordinates  of  the  inner  contour  in  sequence  and  comparing  each 
coordinates’  mapped  outer  coordinate  to  the  last  coordinate  added  to 
the  last  created  continuity  set.  If  that  coordinate  is  within  a 
tolerance  factor  of  the  last  coordinate  added,  it  is  added  to  that  set. 
If  the  coordinate  in  question  is  not  within  tolerance,  a  new  set  is 
created  with  that  coordinate  mapping  as  its  start.  The  tolerance 
factor  used  is  a  ratio  of  the  number  of  coordinates  in  the  outer 
contour  divided  by  the  number  of  coordinates  in  the  inner  contour  times 
a  window  value.  (The  window  value  is  discussed  in  the  next  chapter.) 

'!'o  illustrate  this  continuity  set  assignment,  let  us  refer 
to  the  example  in  Figure  3.11.  Here,  the  tolerance  factor  is  10 
coordinates.  The  last  coordinate  considered  is  inner  coov'dinate  nurnder 
24.  The  next  coordinate  considered  is  coordinate  £5,  which  is  mapped 
to  outer  coordinate  53.  This  coordinate  is  within  the  tolerance  factor 
o'  10  and  is  added  to  the  last  created  continuity  set.  Inner 
coordinate  number  £6  is  mapped  to  outer  coordinate  69.  This  outer 


coordinate  is  outside  of  tolerance  with  the  last  coordinate  added  and 


therefore,  a  new  continuity  set  is  created  with  this  coordinate  mapping 


as  its  start. 

This  initial  step  of  the  continuity  recognition  process  is 
a  fast  method  for  aggregating  coordinate  map  pairs.  In  addition  to 
building  the  initial  continuity  sets  for  the  tentative  mapping,  we  keep 
track  of  the  minimum  and  maximum  mapped  distances  for  each  continuity 
set.  ■^nese  values  are  used  for  merging  continuity  sets  in  the  next 
step  of  the  process. 

The  initial  sets  generated  for  Figures  3.4  and  3. £  are  of 
particular  interest.  This  step  of  the  continuity  procedure  placed  all 
of  the  tentative  mappings  for  the  coordinate  mapping  pairs  for  Figure 
3.4  into  a  single  set.  This  can  be  attributed  once  again  to  the 
contours’  similar  shapes  and  sizes.  On  the  other  hand,  coordinate 
mapping  pairs  for  the  mapping  (1,1)  -  (1,2)  of  Figure  3.6  resulted  in  5 
initial  continuity  sets  with  varying  distance  ranges  (see  Figure  3.12). 

Once  the  initial  continuity  sets  have  been  created  for  a 
contour  pairing,  we  merge  any  sets  that  have  overlapping  mapped 
distance  ranges.  This  merge  process  reduces  the  total  number  of  sets 
and  further  aggregates  the  coordinate  pair  mappings  to  sets  with 
coordinate  number  continuity  and  distance  range  similarity.  In 
re'^erence  to  our  examples,  no  continuity  set  merge  was  required  for 
-igure  3.4  due  to  its  singular  initial  continuity  set.  Figure  3.12 
shows  the  initial  sets  with  distance  ranges  and  the  merged  sets  with 
distance  ranges  for  the  contour  pairing  (1,1)  -  (1,2)  of  Figure  3.6. 
In  that  figure,  the  5  initial  continuity  sets  have  been  merged  into  3 
sets  of  non-overlapping  distance  range. 
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Fig.  3.12  -  Initial  continuity  sets  and  merged  continuity  sets  for 
_ the  contour  pair  (1,1)  -  (1,2)  of  Figure  3.6. 
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Fig.  3.13  -  Bounding  box  overlap  for  exterior  to  exterior  mapping 
Only  the  coordinates  within  the  overlap  area  are  mapped. 
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After  we  have  merged  continuity  sets,  we  need  to  determine 
which  of  those  sets  of  coordinate  mappings  is  the  one  that  should  be 
used  for  connection  formation.  The  choice  is  clearly' the  set  with  the 
smallest  distance  range.  Uith  this  decision,  we  validate  all 
coordinate  pairings  that  are  members  of  this  smallest  distance  set,  and 
cancel  all  other  coordinate  pairings  for  that  set  of  contours, 
b.  Mapping  Cancellation 

The  validated  coordinate  connection  map  for  the  contour 
pair  has  significance  beyond  indicating  which  coordinates  need  to  have 
connection  segments  generated.  It  also  indicates  "filled"  connection 
positions.  By  filled  we  mean  that  once  we  have  formed  connections  to  a 
coordinate  segment  of  a  contour,  that  segment  should  not  be  reused  for 
any  further  mapping  that  occurs  for  the  two  current,  adjacent  planes. 
This  mapping  is  both  checked  and  recorded  at  this  stage  of  the 
algorithm.  Mapping  cancellation  examines  the  coordinate  mappings  for 
which  a  validated  mapping  has  been  assigned.  If  either  of  the  two 
coordinates,  inner  or  outer,  has  been  assigned  to  a  higher  priority 
mapping  for  this  pair  of  planes,  then  that  mapping  is  cancelled.  Once 
these  connections  have  been  struck  from  the  connection  map,  all 
remaining  validated  connections  are  recorded  as  filled. 

An  additional  tasking  of  this  cancellation  process  concerns 
whether  the  mapping  of  either  contour  resulted  in  all  coordinates 
defining  that  contour  being  included  in  the  mapping.  In  that  case,  all 
other  possible  pairings  with  the  completely  mapped  contour  are 
cancelled.  This  is  accomplished  by  zeroing  the  overlap  on  that 
contour’s  row  or  column  of  the  overlap  table. 
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c.  Connection  Formation 

When  the  above  steps  have  been  completed  for  a  pair  of 
contours,  the  remaining  process  of  generating  the'  appropriate  line 
segments  is  relatively  simple.  The  final  coordinate  mapping  for  the 
inner  contour  is  examined  for  continuous  segments  of  validated 
connections.  Mhen  a  continuous  segment  is  defined,  the  beginning  and 
ending  coordinates  of  that  segment  (for  both  the  inner  and  outer 
contours)  are  used  as  boundary  pointers  for  connection  formation.  The 
coordinates  in  between  those  pointers  are  stepped  through  one  at  a  time 
by  a  process  whose  purpose  is  to  generate  the  minimum  area  triangular 
surface  patch,  as  defined  in  our  introduction.  The  surface  patch  is 
formed  by  using  a  line  segment  from  one  contour  as  the  triangle’s  base, 
and  a  coordinate  from  the  other  contour  for  the  triangle’s  third  point. 
The  minimum  area  selection  is  accomplished  by  a  procedure  that  chooses 
the  next  line  segment  between  the  contours  that  is  both  the  shortest 
and  within  the  mapping  specified  for  the  two  contours.  This  is 
identical  to  the  heuristic  used  by  Christiansen  in  CCHRISj.  Differing 
coordinate  rates  between  the  two  contours  are  taken  care  of  by  using 
the  coordinate  ratio  (from  the  continuity  tolerance  factor)  between  the 
contours.  This  ratio  alic'ws  the  process  to  generate  several  line 
segments  emanating  from  a  single  coordinate  when  there  is  a  coordinate 
rata  differential  between  two  mapped  contO'urs.  ~he  lines  generated  by 
this  proicedure  for  Figures  3.4  and  3.  G  are  shown  in  "igures  3.  T  and 
3.7,  respectively. 


Form  the  Coordinats  WaoniriQ:  Exterior 


We  begin  tbe  exterior  to  exterior  mapping  process  at  the  sam 
point  of  the  algorithm  where  we  departed  in  the  description  of  th' 
interior  to  interior  mapping  process.  In  keeping  with  our  order:.' 
criteria  for  mapping  contour  pairs,  we  examine  the  contour  po.i 
requiring  an  exterior  to  exterior  mapping  which  has  the  highest  0'.erla 
percentage  in  the  overlap  table.  fill  remaining  steps  of  the  alg;ri’;h 
are  carried  out  on  this  pair  before  the  next  pair  of  exterio.-'  t 
exterior  contours,  in  largest  to  smallest  overlap  area,  is  o ; ns i det'eo . 

In  Figure  3. 13,  we  are  presented  with  an  enlarged  view  of 
bounding  box  overlap  area  of  the  contour  pairing  (l,i;  -  J'Z,  3.'  ; 
Figure  3.6.  This  area  of  overlap  contains  all  of  the  coiordinates  fro, 
both  contours  which  are  involved  in  the  connection  mapping,  ""he  fi 
ooeration  performed  on  an  exterior  to  exterior  mapped  overlap  pa.-'  i 
the  determination  of  the  set  of  coordinates  in  both  contours  '■hat  i; 
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connect  ion 

the  determination  for  each  coordinate  (in  the  smaller  coor d i -la': e  o-  = 
contour)  of  the  coordinate  in  the  O'ther  co'rtour  coordinate  set  f'at 
the  shortest  distance  away.  This  determinat ion  is  a  simpler  ve'-'sirn  c 
the  distance  minimi cing  process  for  connection  set  assignment 
interior  to  interior  mappings.  The  product  of  this  process  is  f" 

connection  .map  for  the  pair  of  contours.  The  use  of  ccntinuity  sets 
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IV.  ALGORITHM  HEURISTICS  AND  LIMITATIONS 


In  the  pneceeding  chapter,  we  presented  an  explanation  of  our 
algorithm  for  surface  construction.  Particular  attention  was  devoted 
to  our  algorithm’s  handling  of  the  multiple  contours  per  plane  and 
partial  contour  mapping  problems.  It  must  be  emphasized,  however,  that 
our  algorithm  does  not  provide  a  complete  solution  for  all  sets  of 
contour  surface  data.  In  this  chapter,  we  investigate  some  of  the 
limitations  of  our  algorithm.  In  order  to  do  that,  we  must  first 
discuss  the  heuristics  employed  by  that  algorithm. 


A.  HEURISTICS 

Our  algorithm  utilizes  three  heuv'istics  which  are  essential  for  the 
correct  connection  of  planar  contours.  These  heuristics  were  presenteo 
briefly  in  the  last  chapter,  but  we  feel  it  is  necessary  to  explain 
more  fully  their  application  and  interaction  regarding  the  cont  Z'ur 
mapping  problem. 

1 .  Overlap  Percentage  Minimum 

In  step  two  of  our  algorithm,  we  determine  tne  percentage  of 
overlap  between  contours  on  adjacent  planes.  "'’hese  percentages  are 
then  considered  in  a  consistency  check  for  matching  contour  inte"i;'r 
specifications.  The  heuristic  in  question,  the  overlap  percentage 
minimum,  is  applied  in  the  final  phase  of  this  contour  pairing 
procedure.  Contour  pairs  having  an  overlap  percentage  val.e  above  t'-e 
overlap  percentage  minimum,  with  matching  interior  specifications,  are 


designated  for  interior  to  interior  mapping.  Cortoi-u’  pairs  havin: 


non-zero  percentages  below  the  overlap  percentage  on  nimup;,  with 
non-matching  interior  specifications,  are  designated  for  exterior  to 
exterior  mapping,  fill  other  contour  pairs  are  disregarded. 

The  value  we  have  utilized  for  the  overlap  percentage  miriimurn 
is  ten  percent.  We  found,  through  experimentation,  that  the  assignment 
of  this  value  resulted  in  the  greatest  number  of  correct  cont :  .ir 
pairings.  Some  contour  pairs  which  should  be  mapped,  however,  o.;-3 
disregarded  for  mapping  because  of  this  selection  (of  13%)  fo-’  the 
overlap  percentage  minimum.  In  Figure  4.1,  we  are  presented  with  an 
example  of  such  a  situation.  In  that  figure,  we  have  a  pair  o^ 
contours  with  matching  interior  specifications  (HIGH;HIGH),  and  having 
an  overlap  percentage  less  than  ten  percent.  By  our  heuristic,  this 
contour  pair  would  not  be  considered  for  mapping,  and  would  remain 
unconnected. 

One  possible  solution  to  this  problem  would  be  a  vechanisn 
which  used  a  relaxation  procedure  to  force  a  mapping  between  the  pair 
of  contours.  This  mechanism  could  be  selected  by  the  user  to  designate 
contour  pairs  for  mapping  which  would  otherwise  be  disregarded.  1" 
applied  to  the  mapping  situation  of  Figure  4.1,  an  appropriate 
connection  could  be  generated. 

2.  Boundary  Tolerance  Percentage 

The  next  heuristic  to  be  discussed  comes  into  play  in  the 
initial  two  steps  of  our  algorithm.  Specifically,  the  two  operations 
involved  are  the  determination  of  contour  item  two-dimensional  bounding 
box  values,  and  the  usage  of  those  values  for  overlap  determination. 


fis  previously  discussed,  exterior  to  exterior  contour  maooings  are 


Fig.  4.1  -  Example  of  a  contour  pair  which  should  be  mapped,  but 
would  be  disregarded  due  to  overlap  percentage  below  the  minimum. 


indicated  for  pairs  of  contours  with  a  low  percentage  of  overlap  and 


non-matching  interior  specifications.  In  the  initial  developriient  of 
our  algorithm,  we  utilized  the  minimum  and  maximum  X  and  Y  coordinates 
of  the  contour  to  describe  its  bounding  box.  Ue  found,  however,  that 
in  the  majority  of  cases,  these  values  resulted  in  zero  percentage  of 
overlap  between  contours  which  should  be  mapped.  Rn  example  of  this 
limiting  of  bounding  box  values  can  be  seen  in  Figure  4. £,  In  that 
figure,  we  are  presented  with  the  contour  pair  from  -igure  3.  13.  Iri 
this  example,  it  can  be  seen  that  limiting  the  bounding  boxes  for  t^eee 
two  contours  to  their  rescective  minimum  and  maximum  X  and  Y  ccorc  ■.  r  ats 
values  results  in  zero  percentage  of  overlap.  This  is  an 
unsatisfactory  situation  since  the  contours  should  be  mapped. 

To  remedy  this  situation,  we  adjust  the  bounding  box  values  by 
a  percentage  to  promote  mappings  in  situations  similar  to  that 
Figure  4.2.  Once  again,  we  are  presented  with  th=  opportunity  to 

utilize  a  relaxation  procedure,  prompted  via  user  intervent ion,  for 
mapping  situations  not  included  by  this  heuristic.  fl  mechanism  co.ulc 
be  provided  allowing  the  user  to  designate  the  bounding  boxas  fo>' 
individual  contours,  and  thereby  force  a  mapping  between  the  desired 
set  of  contours. 

3.  Tolerance  Mult  ini ier 

In  an  interior  to  interior  mapping  situation,  a  tolerance 
factor  IS  used  for  the  determination  of  the  initial  continuity  set 
assignments.  This  tolerance  factor  is  a  ratio  of  the  number  of 
coordinates  in  the  outer  contou  divided  by  the  number  of  coordinates 
in  the  inner  contour  times  a  window  value.  The  window  value  is  a 
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constant  which  we  found  necessary  for  the  selection  of  appropriate 
mapping  connections.  Ue  chose  to  utilize  a  tolerance  factor  in  this 
step  of  our  algorithm,  as  well  as  in  the  connection  ^orriiation 
procedure,  because  it  provides  an  inexpensive  means  for  restricting  the 
search  space  in  the  selection  of  mapping  connections. 

B.  LIMITOTIONS 

In  the  preceeding  chapter,  we  demonstrated  the  capaoilities  of  :ur 
algorithm,  with  emphasis  on  its  handling  of  the  problems  of  multioie 

9 

contours  per  plane  and  partial  contour  mappings.  We  have  founc, 
however,  that  there  exist  contour  mapping  situations  which  cannot  be 
handled  by  our  algorithm. 

The  first  mapping  situation  concerns  simple  branching  of  one 
contour  on  one  plane  to  two  or  more  contours  on  an  adjacent  plane  (see 
Figure  S.  3) .  In  this  situation,  we  found  that  the  application  o-'  c  ..r 
algorithm  'produces  an  incomplete  contour  mapping  due  to  missing  data. 
One  possible  solution  to  this  mapping  problem  is  the  inclusion  of  a 
crocedure  for  creating  an  introduced  node  similar  to  that  desc^-i'aed  ir 
the  Christiansen  CRef.  £]  paper.  This  special  case  procedure  c;i..ld  te 
selected  automatically,  or  initiated  via  user  interaction. 

The  next  limitation  of  our  algorithm  manifests  itsel^  in  aitaatio's 
where  highly  convoluted  contours,  with  extreme  narrowi ngs,  are  map-igd 
interior  to  interior.  The  problem  here  is  due  to  the  interior  to 
interior  algorithm's  dependence  on  the  overlap  region  bounding  bcx's 
center  coordinate  for  the  tentative  coordinate  mapping.  -or  the 
portion  of  the  contour  near  the  center  coordinate,  the  tentative 


coordinate  is  fairly  good.  For  the  portion  of  a  contour  on  the  other 
side  of  a  narrowing,  where  the  center  coordinate  is  no  longer  central, 
the  tentative  mapping  is  erroneous.  The  problem  ■  comes  when  the 
tentative  mapping  is  so  bad  that  the  continuity  recognition  procedure 
fails,  and  contour  segments  are  incorrectly  left  unconnected. 

The  solution  to  this  problem  is  fairly  simple  and  within  the 
purview  of  our  algorithm.  If  the  convoluted  contour  is  segmented  at 
the  extreme  narrowings,  it  is  possible  to  treat  each  open  segment  of 
the  original  contour  as  a  separate  contour.  Using  the  origirial 
algorithm,  we  can  generate  centers  for  each  new  contour,  and  hence 
coordinate  mappings,  which  result  in  a  more  correct  approx i mat .cr  o" 
the  original  three-dimensional  object.  The  only  capability  lacking 
from  our  present  algorithm  is  a  mechanism  for  partitioning  the  original 
convoluted  contour.  This  mechanism  could  be  either  user  specified  or 
automatic.  The  user  specified  option  is  favored  due  to  the 
computational  expense  involved  for  automatic  contour  segmentation. 

The  next  limitation  also  concerns  interior  to  interioir  contou" 
mappings.  In  situations  where  sections  of  a  contour  tend  to  be  near 
parallel  with  the  vector  drawn  fv'om  the  center  coordinate  of  the  iriner 
contour,  erroneous  mappings  result.  Rn  example  of  this  situation  car 
be  seen  in  Figure  A. 3.  For  those  segments  of  the  outer  contour  whioh 
are  nearly  perpendicular  to  the  tentative  connection  vector,  an 
appropriate  connection  map  is  generated.  fls  the  contour  segrent 
becomes  more  parallel  to  this  vector,  the  tentative  C'onnections 
generated  begin  to  falter. 


erroneous 


tentative 
••  mapping 


Fig.  4.3  -  Example  of  situation  resulting  in  an  erroneous  tentative 
coordinate  mapping  where  contour  segment  becomes  near  parallel 
with  the  tentative  connection  vector. 


Fig.  4.4  -  Example  of  a  situation  where  two  contours  are  mapped 
interior  to  interior  which  would  result  in  an  incomplete  mapping 


The  remedy  to  this  problem  is  very  similar  to  that  for  the  previous 
situation  involving  highly  convoluted  contours  with  extreme  narrowings. 
Segmentation  of  the  original  contour  into  several  open  segments,  which 
could  be  mapped  separately,  would  greatly  improve  the  quality  of  the 
tentative  coordinate  mapping.  Once  again,  user  intervention  is  the 
preferred  method  of  contour  segmentation. 

The  final  problem  situation  to  be  discussed  concerns  interior  to 
interior  mappings  where  the  inner  contour  is  not  contained  in  the  o.iter 
contour.  This  situation  would  result  from  contour  data  taken  from  a 
torus,  such  as  a  doughnut.  fln  example  is  illustrated  in  Pigoire  4.4. 
The  problem  with  this  mapping  situation  results  from  the  use  of  the 
tentative  connection  vector  emanating  from  the  center  of  the  iomer 
contour.  Since  the  center  coordinate  of  the  inner  contour  is  displaced 

from  the  center  coordinate  of  the  outer  contour,  tentative  mappings  a-2 
generated  only  for  that  section  of  the  outer  contour  which  is  the 
same  side  of  the  tentative  connection  vector  (see  Figure  4.4).  '"he  net 
result  is  a  partial  mapping  of  two  contours  which  should  be  totally 
connected. 

fl  practical  solution  to  this  mapping  problem,  which  could  fcs 
readily  adapted  to  our  algorithm,  is  described  in  the  Christ iarsen 
paper  CRef.  £3.  In  mapping  situations  wnere  contours  to  be  mapped  are 
not  mutually  centered,  Christiansen  recommends  a  translation  procec‘..re 
onto  a  unit  square,  centered  at  (0,0).  The  principle  of  this  process 
is  to  translate  the  two  contours  in  such  a  manner  that  they  become 
mutually  centered  within  the  unit  square.  Ppplication  of  the  interior 
to  interior  algorithm  at  this  ooint  would  result  in  the  desirsc 


mappings.  Tentative  mappings  would  be  generated  for  the  contouv's’ 
original  coordinates,  thus  allowing  the  appropriate  connections  to  oe 
formed  in  the  final  step  of  the  algorithm. 

C.  SUMMARY 

It  has  been  the  ourpose  of  this  chapter  to  investigate  the 
limitations  of  our  algorithm,  and  provide  practical  solutions  where 
possible.  Additionally,  to  a  lesser  degree,  the  heuristics  employed  by 
our  algorithm  have  been  explained  to  increase  the  understand ing  o'®  t^e 
reader.  We  feel  that  our  algorithm  provides  a  solution  to  the  multiple 
contours  per  plane  and  partial  mapping  problems,  but  must  concede  that 
it  is  not  a  total  solution  to  the  problem  of  surface  construction  -^rcm 


planar  contour  data. 


APPENDIX 


PSEUDO-CODE  ALGORITHM  DESCRIPTION 


A  pseudo-code  description  of  our  surface  construction  algorithni 
follows: 

FACE_IT 

Input  the  coordinates  for  two  adjacent  planes.  Make  a  local  copy 
of  the  coordinates. 

DELINEATE_INVENTORY 

< 

Take  inventory  of  the  contours  in  the  coordinate  sets.  This 
inventory  determines  the  total  number  of  contours  for  each  plane 
and  records  where  each  contour  begins  and  ends. 

> 

TYPE_INVENTORY 

We  determine  the  contour  type  of  each  contour  in  each  plane. 
There  are  three  possible  types:  CLOSED_LCCP,  OPEN_SEGMENT,  and 
SINGLE_POINT. 

> 

BOUND_INVENTQRY 

Determine  the  rectangular,  two-dimensional  boundary  of  each 
contour.  Increase  those  boundaries  by  a  constant  to  increase 
the  possibility  of  detecting  appropriate  exterior  to  exterior 
mappings. 

> 

I NTER I □R_DETERM I NAT I ON 
■C 

Determine  whether  the  interior  of  each  contour  is  HIGH  or  LOW 
valued  with  respect  to  the  current  contour  level.  This  valus 
can  be  assigned  interactively  in  cases  where  the  information 
to  make  this  determination  is  not  available.  These  values 
are  used  in  a  consistency  check  for  selection  of  contour  pairs 
for  mapping. 

> 

OVERLAP_DETERM I NATION 
{ 

Compute  the  overlap  table  for  the  contours  of  both  planes.  The 
values  in  the  table  are  the  percentage  of  overlap  for  each 
possible  contour  pair  on  the  adjacent  planes.  If  there  is  no 
overlap,  a  value  of  0.0  is  recorded. 


h 


Contour  mapping  types  are  also  assigned  at  this  step  of  the 
algorithm.  Contour  pairs  with  a  HIGH  percentage  of  overlap, 
matching  interior  specifications  (HIGH;HIGH,  LOW:LOW,  or 
INDETERMINATE: anything)  are  assigned  interior  to  interior  type 
mapping.  Those  pairs  with  a  non-zero  overlap  percentage,  belcw 
10%,  with  non-matching  interiors  are  assigned  exterior  to 
exterior  mappings.  All  other  contour  pairings  are  zeroed. 


CONNECT I ON_DETERM I NAT I ON 

< 

This  step  of  the  algorithm  orders  the  pairs  to  be  rnaoped,  and 
forms  connections  for  the  assigned  types  of  contour  mappings. 
This  step  is  detailed  below. 

> 


>  /■*  end  of  FACEIT  #/ 


CONNECT I ON_DETERM I NAT I ON 
{ 

while  .true. 

{ 

Find  the  largest  overlap  percentage  in  the  overlap  table.  If  the 
largest  value  =0.0  then  QUIT. 

If  the  contour  mapping  indicated  by  this  largest  overlap  value  is 
exterior  to  exterior  .... 

EXTERIOR_TO_EXTERIOR_MAPPING 

< 

Determine  the  set  of  coordinates  in  each  contour  that  are  in  the 
overlap  area. 

For  the  contour  of  the  overlap  pair  that  has  the  least  number  of 
coordinates,  find  the  minimum  distanced  coordinate  of  the  other 
contour. 

Assign  all  coordinates  within  the  overlap  region  to  the  connec¬ 
tion  set. 

>  /*  end  of  EXTERIOR_TO_EXTERIOR_MAPPING  */ 

else 

/*  perform  an  interior  to  interior  mapping  ■*/ 

INTERIQR_TO  INTERIOR  MAPPING 


Determine  which  contour  of  the  pair  is  interior.  This  assignmen 


IS  based  upon  which  contours’  bounding  box  is  smallest. 


Compute  the  center  coordinate  of  the  inner  contour’s  bounding 
box.  Check  to  make  sure  that  this  point  is  inside  the  contour. 
If  it  is  not,  the  contour  needs  to  be  partitioned". 

For  each  coordinate  of  the  inner  contour,  determine  the  coord¬ 
inate  of  the  outer  contour  which  is  closest  to  a  vector  drawn 
from  the  center  coordinate  through  the  coordinate  of  the  inner 
contour.  Store  the  coordinate  as  the  connection  map  coordinate 
for  the  inner  contour.  Plso,  record  the  mapped  distance  from 
each  inner  coordinate  to  its  mapped  outer  coordinate. 


RECQGN I Z  E_CONT I NU I TY 

/*  Determine  continuity  sets  in  the  two  contours  using  the 
the  connection  map  and  associated  distances.  */ 


I N I T I AL_CONT I NU I TY_SETS 
■C 

Assign  the  coordinates  of  the  connection  map  to  a  con¬ 
tinuity  set  based  upon  whether  each  consecutive  coordinate 
is  within  a  coordinate  tolerance  factor.  ’"nis  tolerance 
factor  is  a  ratio  of  the  number  of  coordinates  in  the 
outer  contour  divided  by  the  number  of  coordinates  in  the 
inner  contour. 


I N I T I AL_SET_D I STANCE.RANGES 

< 

Determine  the  minimum  and  maximum  distance  ranges  for  each 
of  the  continuity  sets. 

> 

CONT I NU I TY_SET_MERGE 

< 

Merge  any  continuity  sets  that  have  overlapping  distance 
ranges,  maintaining  the  distance  range  for  any  merged  set. 

> 

CONNECT I ON_SET_ASS I GNMENT 

Assign  coordinate  connections  for  the  coordinates  of  the 
merged  continuity  set  that  contains  the  smallest  distance. 
All  other  continuity  sets  are  left  unconnected. 

> 

>  /*  end  of  RECOGNIZE_CONTINUITY  */ 

/*  end  of  INTERIOR  TO  INTERIOR  MAPPING  */ 


h*-  **-  .“'i  ,*•  w"*  »"•  V**  •.**•*'»* 


MflPP I NG.CflNCELLflT I ON 

< 

Examine  the  coordinate  maopings  for  which  a  connection  has  been 
assigned.  If  either  of  the  two  coordinates,  inner  contour  or 
outer  contour,  has  been  used  in  a  previous,  higher  priority 
mapping  for  this  pair  of  planes,  that  coordinate  mapping  is 
cancelled.  Once  these  filled  connections  have  been  struck  from 
the  connection  map,  all  remaining  validated  connections  are 
recorded  as  filled. 


CONNECT! DN_FORMflT I ON 

< 

Generate  the  connections  for  the  validated  coordinate  map.  "^hi 
is  accomplished  by  stepping  through  the  connection  map  ard 
forming  coordinate  connections  where  indicated.  In  between 
coordinates,  those  not  directly  mapped  but  within  the  tolerance 
factor  for  the  connection  mapping,  are  also  added  to  the  oictur 
The  goal  of  the  connection  process  is  to  form  minimum  arsa’d 
triangular  surface  patches. 

> 

>  /•*  end  while  .true,  of  CONNECTION  DETERMINPTION  */ 


>  /♦  end  of  CONNECTION  DETERMINATION  »/ 
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